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Abstract 

We present a new method to investigate a class of diagrams which generalizes the 
sunset topology to any number of massive internal lines. Our attention is focused on 
the computation of the spectral density of these diagrams which is related to many- 
body phase space in D dimensional space-time. The spectral density is determined 
by the inverse i^-transform of the product of propagators in configuration space. The 
inverse if-transform reduces to the inverse Laplace transform in any odd number of 
space-time dimensions for which we present an explicit analytical result. 
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(a) (b) 

Figure 1: (a) a sunset diagram with three different masses nii, m2 and 
(b) general topology of the class of water melon diagrams 



Recently there has been renewed interest in the calculation of sunset-type diagrams with 
massive internal lines (as some of the latest refs. see for example [|l], 0, ^, ^). The simplicity 
of the sunset-type topology has a special appeal since sunset-type diagrams can be regarded 
as a laboratory for testing new methods in the computation of diagrams with massive 
particles (see e.g. p). For phenomenological applications one needs to avail of efficient 
numerical algorithms to compute the sunset-type diagrams (see e.g. Q). In ||^ we have 
derived a one-dimensional integral representation for a class of diagrams that generalizes 
the sunset topology to any number of internal lines (we call them water melon diagrams). 
We believe that the method described in [0 and in the present paper completely solves the 
problem of calculating the class of sunset-type diagrams. 

In this note we present a new representation for the spectral density of water melon 
diagrams through an integral transformation in configuration space. It is given by an 
inverse i^'-transform of the product of propagators of internal lines that constitutes the 
water melon diagram. The representation simplifies to the inverse of the standard Laplace 
transform in any odd number of space-time dimensions for which we present a concise 
analytical result. The spectral density of water melon diagrams is related to the many- 
body phase space integral in D dimensional space-time for which we find novel results. 

The sunset diagram proper represents the leading order perturbation theory correction 
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to the propagator in (^"^-theory, i.e. it is a two-loop insertion with three internal lines (see 
Fig. 1(a)). A straightforward generalization of this topology (water melon diagram) is a 
correction to the propagator in the ^'^"'"^-theory that contains n loops and n + 1 internal 
lines (see Fig. 1(b)). As an example take the leading quantum corrections of higher order 
chiral perturbation theory for pseudoscalar mesons where such diagrams occur. Thus we 
describe as a generic case of the water melon diagram the correlator of two field monomials 
jn(yX) =: (/>i,^i . . . 0n,/x„: taken in the normal ordering form where is a derivative of the 
field (pi with a multi index /i = {/ii, . . . ,/ip}. We leave aside trivial generalizations that 
would result from the inclusion of derivatives and the emission/absorption of additional 
external particles with vanishing momentum. These would merely result in the doubling 
of the propagators to which they are attached and thus are easily included in our general 
approach. Returning to the scalar case the general representation of the scalar polarization 
function corresponding to water melon diagrams is given by 

n 

Il{x) = l[D{x,m,) (1) 

i=l 

where D{x,m) is the configuration space propagator of a particle with mass m in D- 
dimensional (Euclidean) space-time ^ 

" (2^ J p^pi^ + m2 " (2^)A+ia;2A ^ ) 

where D = 2X + 2, and Kx{z) is a McDonald function (a modified Bessel function of the 
third kind, see e.g. |^). When we refer to a specific space-time dimension D later on 
we add the space-time index to the propagator function, i.e. we write Do{x,m) for the 
propagator function. The propagator depends only on the length x of the D-dimensional 
vector (x^), x = ^/x^x^ (a times we shall also use the notation x for the vector itself). 

The representation Eq. (|l]) contains the complete information about the class of water 
melon diagrams and is the starting point for all further investigations. 

For some applications one needs the Fourier transform of the polarization function 



which is defined by 

U{p) = J U{x)e'P-^'d''x = I {TUx)j^me'P-^U''x. (3) 

The analytic structure of the Fourier transform 11 (p) is determined by the Kallen-Lehmann 
representation 



n(p) = / ' ^ - (subtractions) (4) 



p{s)ds 
h s + p"^ 

with the appropriate number of subtractions for the regularization of ultraviolet diver- 
gences, if necessary. The spectral density p(s) or the discontinuity across the cut in the 
complex plane of the analytic function Il{p) is of interest for various physics applications. 
We mention, though, that in the Euclidean region the function Il{p) itself also has some 
applications. For these applications the values of Il{p) can easily be found from the rep- 
resentation Eq. or directly from Eq. (§) 0]. Here we shall focus our attention on the 
spectral decomposition of the water melon diagrams which in turn is connected with the 
particle content of the underlying theory. 

The analytic structure of the correlator Il{x) (or the spectral density of the correspond- 
ing polarization operator) can be determined directly in configuration space without having 
to compute its Fourier transform first. The dispersion representation (or the spectral de- 
composition) of the polarization function in configuration space has the form 

POD 

U{x) = / p{s)D{x,s)ds (5) 
Jo 

where D{x, s) is the scalar propagator with the mass m = ^/s. This representation was 
used for sum rules applications in [0, |ll| where the spectral density for the two-loop 



sunset diagram was found in two-dimensional space-time |]T2l. With the explicit form 



of the propagator in configuration space Eq. (|^) the representation in Eq. (^ turns into 
a particular example of the Hankel transform, namely the iC-transform [IB, 13 1. Up to 
inessential factors of x and ^/s, Eq. (Kf) reduces to the generic form of the X-transform for 
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a conjugate pair of functions / and g, 

roo 

9{y) = / f{x)K^{xy)y/xydx. (6) 
Jo 

The inverse of this transform is known to be given by 

]^ pc+ioo 

f{x) = — g{y)I^{xy)y/xydy (7) 

m Jc-ioo 

where Iu{x) is a modified Bessel function of the first kind and the integration runs along 
a vertical contour in the complex plane to the right of the right-most singularity of the 
function g{y). In order to obtain a representation for the spectral density p{s) of water 
melon diagrams in general D-dimensional space-time one needs to apply the inverse K- 
transform to the particular case Eq. (^. One has 

(2,7t)'^ fc+ioo 

m^p{m^) = ^—^ / Il{x)x^^^Iximx)dx. (8) 

I Jc-ioo 

The transform Eq. (Q) solves the problem of determinating the spectral density of water 
melon diagrams by reducing it to the computation of a one-dimensional integral for a 
general water melon diagram with any number of internal lines and different masses. Below 
we discuss some technicalities concerning the efficient evaluation of the contour integral. 

The analytic structure of the correlator in Eq. (|l]) is now explicit using the represen- 
tation given by Eq. (||) and exhibits the distribution nature of the spectral density p(s). 
For odd-dimensional space-time the representation in Eq. (|^) reduces to the Laplace trans- 
formation. To obtain the spectral density (the function f{x) in this particular example) 
one can use Eq. (|^). For energies below threshold it is possible to close the contour of 
integration to the right. With the appropiate choice of the constant c as specified above 
the closed contour integration gives zero due to the absence of singularities in the relevant 
domain of the right semi-plane. By closing the contour of integration to the left and keep- 
ing only that part of the function I^iz) which is exponentially falling for Re(2;) < one can 
obtain another convenient integral representation for the spectral density when the energy 
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is above threshold. The only singularities within the closed contour are then poles at the 
origin (in odd-dimensional space-time) and the evaluation of the integral can be done by 
determining the corresponding residues. These are purely algebraic manipulations, the 
simplicity of which also explain the simplicity of the computations in odd-dimensional 
space time. For a small number of internal lines n the spectral density can also be found 
by using the convolution formulae for the spectral densities of a smaller number of particles 



(see e.g. [1T5|). For large n the computations described in [|15l become quite cumbersome 
and the technique suggested in the present paper is much more convenient. For example, 
in odd-dimensional space-time our techniques give the result for water melon diagrams in 
terms of elementary functions. In three space-time dimensions our results can be used to 
compute the phase-space integrals for particles in jets where the momentum along the di- 
rection of the jet is fixed |T^. Another application can be found in three dimensional QCD 



which emerges as the high temperature limit of the ordinary theory of strong interactions 
for the quark-gluon plasma (see e.g. [jl^ |18|, 

For even-dimensional space-time the analytic structure of 11 (x) in Eq. (|l]) is more com- 
plicated: there is a cut along the negative semi-axis which prevents a straightforward 
evaluation by simply closing the contour of integration to the left. The discontinuity along 
the cut is, however, well known and includes only Bessel functions that appear in the 
product of propagators for the polarization function. Therefore the representation (P) is 
essentially equivalent to the direct analytic continuation of the Fourier transform but 
may be more convenient for numerical treatment because there is no oscillating integrand 
in®. 

Note also that the expression given by Eq. (|l|) can have non-integrable singularities 
at small x for a sufficiently large number of propagators when D > 2 |^. Therefore 



the computation of its Fourier transform requires regularization or subtractions [0. The 
spectral density itself is finite (the structure of the water melon diagrams is very simple 
and there are no sub divergences when one employs a properly defined i?-operation ||21|| ) 
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and thus requires no regularization. The exphcit representation for the spectral density 
within the more traditional approach reads M 



2-Ki Discn(p)|,n2 



1 4 V 2< ; V , 



d^, (9) 



where one has taken the discontinuity of the Fourier transform across the physical cut for 
= —s ± iO, s > 0. The analytic continuation is done according to [§] 

Kx{z) = ^ei^^i7«(zz) (10) 

with H^''^\z) being the Hankel functions, H^\z) = {H^^\z))* for real z and A. This is an 
alternative representation of the spectral density, and in some instances this representation 
can be more convenient for numerical treatment. 

In the following we present some explicit examples of applying the technique of comput- 
ing the spectral density of water melon diagrams on the basis of integral transformations 
in configuration space. 

First a remark about the mass degenerate one-loop case is in order. All necessary 
integrals (both for the direct and the inverse i^'-transform) involve no more than the 
product of three Bessel functions which can be found in a standard collection of formulas 
for special functions (see e.g. |^). The spectral density in D-dimensional space-time (for 
two internal lines with equal masses) can be computed to be 

pis) = ^ — — ; > s>Am\ (11) 

^ 24A+%^+lr(A + i)v/i ^ ' 

This formula is useful since it can be used to test the limiting cases of more general 
results. Similarly one can derive a general D-dimensional analytical result for the one-loop 
nondegenerate mass case. 
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For three-dimensional space-time we have the dispersion representation 

Il(x) = / p{m^)D^ix^'m)dm'^ = / pim^)— dw?' (12) 

JO Jo 4:7TX 



with the three-dimensional scalar propagator 



y/mxKi{mx) ^-mx 



One can invert Eq. ( |T^ ) and obtains 

]^ rc+ioo 



I PC+'LOO 

2mp{m^) = — : / 47rxn(x)e"^c/x (14) 
which is a special case of Eqs. and (||) with 

^iW = \/^sinh(z) (15) 

where one only needs to retain the piece in the hyperbolic sine function. The solution 
given by Eq. (jM]) has the appropriate support as a distribution or, equivalently, as an in- 
verse Laplace transform. It vanishes for m < M = I]r=i since the contour of integration 
can be closed to the right where there are no singularities of the integrand. Recall that 
for large x with Re(a;) > the asymptotic behaviour of the polarization function Il{x) is 
governed by the sum of the masses of the propagators and reads 

n(a;) ~ exp(-xM). (16) 

For m > M one can close the contour to the left and then the only singularities of 11 (x) 
are the poles at the origin of 11 (x) since it is a product of the propagators of the form of 
Eqs. (|l]) and (p^ . The integration in Eq. ( |I4D then reduces to finding the residues. Indeed, 

n -Mx 

U{x) = l[D,{x,m,) = -——, (17) 
~=i [Airxp 

and Eq. (0) gives 

I i-c+ioo \ j-c+ioo Q{m-M)x 



2mp(m^) = / 47rxn(x)e"'^'dx = r--— / 

' 27TiJc-ioo ^ ^ 27rz 47r "-iic 



27rz(47r)"-i Jc-ioc 
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-dx. 



At m > M one closes the contour of integration to the left and computes the residue at 
the origin to obtain (n > 1) 

im — M)""^ 

2nipim^) = , -9{m - M). (19) 

'^^ ' (47r)"-i(n - 2)! ^ ^ ^ ^ 

This also explains the simplicity of the structure of the spectral density in odd numbers of 
dimensions of space-time when traditional means are used [|^. In five- dimensional space- 
time Eq. (|1^) is applicable almost without any change because the propagator now reads 

3 

{mx)'^ Ki{mx) ^-mx 
D,{x,m)= .g/, 3 =-^il + mx). (20) 

Compared to the three-dimensional case the only additional complication is that the degree 
of the order of the pole at the origin is changed and that one now has a linear combination 
of terms instead of the simple monomial in three dimensions. 

In even number of dimensions one is dealing with a genuine K-transform. We discuss 
in some detail the important case of four- dimensional space-time. For D = 4 (A = 1), 
Eqs. (@) and (||) give 

n(a;) = / pim')D^{x,m)dm^ = f p{m^)'^^^^^^^$^dm\ (21) 
J J Atx^x^ 

and Eq. (|) is written as 

2mp{m^) = — / AiT^x'^U{x)Ii{mx)dx. (22) 

All remarks about the behaviour at large x apply here as well. However, the structure 
of singularities is more complicated than in the odd-dimensional case. In addition to the 
poles at the origin there is a cut along the negative semi-axis that renders the computation 
of the spectral density more involved. The cut arises from the presence of the functions 
Kiiniix) in the polarization function n(x). Also the asymptotic behaviour of the func- 
tion Ii{z) is more complicated than that of Ii/2{z). In particular the extraction of the 
exponentially falling component on the negative real axis is rather tricky. Incidentally, the 
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fall-off behaviour of the function /i (z) on the negative real axis can be taken as an illus- 
tration of Stokes' phenomenon of asymptotic expansions (see e.g. 0). While the analytic 
structure of the representation is quite transparent and the integration can be performed 
along a contour in the complex plane, there are some subtleties when one wants to obtain 
a convenient form for numerical treatment 0. 

After closing the contour to the left (for m > M) using the appropriate part of the 
function Ii{z) we obtain 

pc+ioo 

m / x'^I\.{x)Ii{mx)dx 

J c—ioo 

Ki{mr)r^ (n(e^^r) + n(e"^''r)) dr + 2 ^ Ki{mryU{r)dr (23) 

+ / {i7rli{mz) + Ki{mz))Il{z)z'^dz + / {iTTli{mz) — Ki{mz))Il{z)z'^dz 
JC- Jc+ 

where the contours C+ and C_ are semi-circles of radius e around the origin in the upper and 
lower complex semi-plane, respectively. For practical evaluations of n(e^'"^r) the following 
rule for the analytic continuation of McDonald functions is used: 

Ki{e^^^mr) = —Ki{mr) ^ iTrIi{mr), mr > 0. (24) 

Let us add a few remarks on the final representation Eq. (|23|) which is in a suitable 
form for numerical integration. We have introduced an auxiliary regularization in terms 
of a circle of finite radius e which runs around the origin with its pole-type singularities. 
The spectral density is independent of e, and the parameter e completely cancels in the full 
expression for the spectral density as given by Eq. (pSf). This must be so since the spectral 
density is finite for the water melon class of diagrams. Eq. (|23D contains no oscillatory 
integrands (cf. Eq. @), and the integration can be safely done numerically. Thus Eq. ( ^3]) 
is a useful alternative representation for the spectral density. In practice the integration 
over the semi-circles is done by expanding the integrand in z for small z and keeping 
only terms singular in e. The expansion requires only a finite number of terms and is a 
purely algebraic operation. Then the singularity in e exactly cancels against those of the 
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remaining integrals. This cancellation can also be done analytically leaving well defined 
and smooth integrands for further numerical treatment. 

Next we turn to the threshold behaviour of the spectral density. Using the results of 
the above analysis one finds 



where Z'i(mjr) is either /i(mjr) or Ki{mir). The convergence at large r > is controlled 
by the factor exp (— (m — M)r) as in Eq. (|l^), and the corresponding expansions in the 
variable m — M in the region m ^ M can be easily constructed. 

To summarize, we have described a novel technique to compute the spectral density of 
any water melon diagram which reduces the calculation to a one- dimensional integral with 
well known functions in the integrand. Any tensor structure and form factor structure can 
be easily added without modification of the basic formulae. Different regimes of behavior 
with respect to mass/momentum expansions can easily be written down in terms of only 
two dimensionful parameters (one mass and one momentum or two masses) while other 
parameters are taken as small. The threshold behaviour of the spectral density can be easily 
investigated based on such a representation. Explicit analytical formulae (with even that 
last integration being performed) are found in odd space-time dimensions. This allows one 
to compute any n-particle phase space integral for any kind of particles. The only reason 
that prevents the final integration to be explicitly done in even number of dimensions is that 
products of Bessel functions are encountered which cannot be integrated in closed form. 
This forces one to use numerical integrations in the even-dimensional case except for very 
special cases. Even then the analytic structure of the solution is completely determined 
which allows for an efficient and reliable numerical treatment of the problem. 

All in all we conclude that the problem of computing the class of multi-loop massive 
water melon diagrams can be considered to be completely solved. 




(25) 
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